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Abstract 

The Standard Model predictions for the hadronic and leptonic electric dipole moments 
(EDMs) are well far from the present experimental resolutions, thus, the EDMs represent 
very clean probes of New Physics effects. Especially, within supersymmetric frameworks 
with flavor-violating soft terms large and potentially visible effects to the EDMs are 
typically expected. In this work, we systematically evaluate the predictions for the EDMs 
at the beyond-leading-order (BLO). In fact, we show that BLO contributions to the 
EDMs dominate over the leading-order (LO) effects in large regions of the supersymmetric 
parameter space. Hence, their inclusion in the evaluation of the EDMs is unavoidable. 
As an example, we show the relevance of BLO effects to the EDMs for a SUSY SU (5) 
model with right-handed neutrinos. 



1 Introduction 

The Standard Model (SM) of elementary particles has successfully passed all the elec- 
troweak tests at the LEP and also all the low-energy flavor physics tests. On the other 
hand, there is a general agreement that the SM has to be considered as an effective field 
theory, valid up to some unknown scale of New Physics (NP). A natural solution of the 
hierarchy problem would point towards a scale of new physics close to the TeV, an energy 
scale that will be explored by the LHC 
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Besides the direct NP search at the coUiders (the so-called high- energy frontier), a com- 
plementary tool to shed light on NP is provided by high-precision low-energy experiments 
(the so-called high-intensity frontier) specially to determine the symmetry properties of 
the underlying NP theory. The hadronic uncertainties and the overall good agreement of 
flavor-changing neutral current (FCNC) data with the SM predictions prevent any con- 
clusive evidence of NP effects in low-energy precision tests of the quark sector. However, 
a closer look at several CP- violating observables indicates that the CKM phase might not 
be sufficient to describe simultaneously CP violation in K, Bfi and Bg meson decays. In 
particular: 

• The values of sin 2/5 extracted from Penguin-dominated modes, i. e. 3^ 4>Ks, f]'Ks-, 
tt'^Ks, ^Ksi KsKsKs etc., are significantly lower than the corresponding value from 
the tree-level process Ed 'ipKs, i.e. (sin 2/3)^;^^ = 0.680 ± 0.025 [1]. 

• CP violation in the Bd — Bd system, i.e. (sin 2/9),/,k<,, combined with the AMd/AMs 
constraint, appears insufficient to describe the experimental value of ex within the 
SM, after the new SM value for €k is taken into account [2] . Similarly, a simultaneous 
description of €k and AMd/AM^ within the SM implies sin 2/3 = 0.88 ± 0.11 [3] [2], 
well far from the measured (sin 2/3)^^-^. 

• There are some hints for the asymmetry in Bg —>■ %l)(f), S^^, to be significantly larger 
than the SM value S^^ 0.04 [4]. 

• There are also other interesting tensions observed in the data, as the rather large 
difference in the direct CP asymmetries Acp {B- K-n°) and Acp(Rd K-n+) 
and certain puzzles in Bd — > ttK decays. 

All the above tensions might be accommodated through the introduction of new sources 
of CP violation in a given NP scenario. 

We know that some sources of CP violation, in addition to the unique CKM phase, are 
required to reproduce the observed baryon asymmetry of the universe. In this context, 
it seems of particular interest the study of those low-energy observables that are quite 
sensitive to new sources of CP violation, as the Electric Dipole Moments (EDMs). In fact, 
the EDMs offer a unique possibility to shed light in NP, given their strong suppression 
within the SM and their high sensitivities to NP effects. 

The minimal supersymmetric SM (MSSM), that is probably the most motivated model 
beyond the SM, exhibits plenty of CP-violating phases [5] able to generate the EDMs at 
an experimentally visible level [6-8]. After SUSY-breaking terms are introduced, new 
CP- violating sources may naturally appear through i) flavor-conserving terms such as the 
B^ parameter in the Higgs potential or the A terms for trilinear scalar couplings and ii) 
flavor- violating terms such as the squark and slepton mass terms. It seems quite likely 
that the two categories i) and ii) of CP violation are controlled by different physical 
mechanisms, thus, they should be distinguished and discussed independently. 

In the case i), non- vanishing quark EDMs and lepton EDMs are generated already 
at one-loop level. It is well-known that the flavor-conserving phases 0a and 0^ have 
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to be very small, at the level of O(10~^^~^^), if a low-energy SUSY is realized around 
the electroweak scale. Obviously, some mechanisms are required to suppress them in a 
natural way. The SUSY CP problem can be avoided in a number of SUSY scenarios. 
The gaugino/gauge mediation models for SUSY breaking arc candidates when the A and 
B/j, parameters are derived by the renormalization group (RG) effects. We also remind 
the "effective SUSY scenarios" [9], the "Split SUSY hmit" [10], or the "Dirac gaugino 
models" [11]. 

If the flavor- conserving terms do not have CP-phases, CP might be violated only by 
those soft terms that additionally change the flavors, i.e. the case ii). The flavored 
EDMs, induced by the flavor- violation, turn out to be naturally suppressed because of 
the smallness of the mixing angles regulating the flavor transitions. In such a case, the 
EDMs and flavor-changing processes are intimately related as they both come from the 
same sources and their correlated study can represent a useful tool to test or to falsify 
this scenario. 

In this paper, we analyze in detail the predictions for the flavored EDMs and Chromo- 
EDMs (CEDM), at the beyond-leading-order (BLO) level within SUSY theories. As 
already stressed in Refs. [12,13], the inclusion of BLO effects is essential for a correct 
evaluation of the flavored EDMs. In fact, BLO effects do not just represent a correction 
to the leading-order (LO) contributions, but, they can dominate over the LO ones in a 
broad region of the SUSY parameter space. 

At the LO, the hadronic EDMs are generated by the one- loop exchange of gluinos 
(g) and charginos {x"^) with squarks. The dominant BLO contributions are computed by 
including all the one-loop induced (tan /9-enhanced) non-holomorphic corrections for the 
charged Higgs (i?^) couplings with fcrmions and for the x'^/g couplings with fermions- 
sfermions. The above effective couplings lead to the generation of effects to the 
(C)EDMs, absent at the LO, via the one-loop if^/top-quark exchange. Moreover, the 
chargino contributions, suppressed at the LO by the light quark masses, are strongly 
enhanced at the BLO by the heaviest-quark Yukawa couplings. Finally, also the gluino 
effects receive large BLO contributions that are competitive, in many cases, with the 
LO ones. As a result, BLO effects dominate over the LO effects in large regions of the 
SUSY parameters space and their inclusion in the evaluation of the flavored (C)EDMs is 
mandatory. 

The organization of the paper is as follows. In Section 2, we review the LO contri- 
butions to the flavored EDMs. Here, we also discuss the Jarlskog invariants, to which 
the (C)EDMs are proportional to. The invariants are a useful tool to classify the CP 
violation in the flavor-changing SUSY-breaking terms. In Section 3, we show the BLO 
contributions to the flavored EDMs and we compare them with the LO contributions. 
For completeness, the BLO contributions to the electron EDM and the CP-violating four- 
Fermi operators are also given in the section. In Section 4, the EDMs in the SUSY GUTs 
are discussed. Section 5 is devoted to the conclusions. 

In Appendix A, we show the derivation of various effective couplings induced by the 
tan 5-cnhanced radiative corrections with explicit CP phase dependence. The effective 
vertices in the mass insertion approximation is shown there. The loop factors used in 
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Appendix A are given in Appendix B. In Appendix C we present the effective vertices in 
the mass eigenstate basis necessary for the numerical calculation. Appendix D contains 
the formulae for the (C)EDMs in the mass eigenstate basis, using the effective vertices. 
In Appendix E and F the loop functions used in this paper are defined. 



2 Jarlskog Invariants and Flavored EDMs at the LO 

The effective CP-odd Lagrangian, which contributes to the (C)EDMs, is given as 



i=u,d,s,e,fi i=u,d,s 

r'^G^G'^^^'G^^'^ + ^ a, (^,V^,)(V^,275V^,) + ■ • • , (1) 

where F^y and are the electromagnetic and the SU{'i)c gauge field strengths, respec- 
tively. The first term of Eq. ([1]) is the well-known QCD theta term. The second and third 
terms of Eq. ([T]) are the fermion EDMs and CEDMs, respectively, while the second line 
of Eq. ([T]) contains dimension-six CP-odd operators, such as the Weinberg operator and 
the CP-odd four-Fermi interactions. 

The QCD theta parameter is tightly constrained by the neutron EDM at the level of 
6 < 10"(^~^''); this naturalness problem is commonly referred to as the strong CP problem. 
One natural way to achieve the required suppression for 9 is to impose the Peccei-Quinn 
symmetry, since the axion field makes 9 dynamically vanishing. Under the above assump- 
tions, the quark (C)EDMs are very suppressed in the SM as they are generated only at 
the three-loop level by the phase of the CKM matrix. Long-range effects to the neutron 
EDM, arising at the two-loop level, are still far below the current bound. In this paper 
we assume the Peccei-Quinn symmetry, for simplicity. 

Among the various atomic and hadronic EDMs, a particularly important role is played 
by the thallium EDM (liTi) and by the neutron EDM [dn) that can be estimated as [14-16] 

rf^i = -585 - e 43 GeV Cf , (2) 

where Gf^ has been evaluated as [17] 

(0) _ 29MeV K X 220 MeV 66 MeV(l - 0.25/t) 

^ S ~ ^de 1" ^se 1" ^be ; ['->) 

rrid rus m^, 

with ft ~ 0.5 [18]. Moreover, dn can be estimated as [19,20] 

dn = {l± 0.5) [l.4 {dd - 0.25 du) + 1.1 e (rf;^ + 0.5 <)] . (4) 

In our analysis, we will use the above formulae for the evaluation of the (C)EDMs0. 



^ It has been argued in Ref. [21] that the evaluation of the neutron EDM stiU suffers from sizable 
uncertainties even when evaluated within the QCD sum rules approach. 
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The strong suppression of the quark (C)EDMs in the SM could be understood by 
considering Jarlskog invariants (JIs) which represent a basis-independent measure of CP 
violation. In the SM, we have the following JIs, 

Js'i, = im{y,y„[y,,y„]y,},,, 

Js^ = im{y„y„[y,,y„]y,}n (5) 

for down- and up-typc quarks, respectively. Here, {yd{u))ij is the down(up)-type quark 
Yukawa couplings and Yd(u) = y\{u)y<i{u)- Since the EDMs and CEDMs of the z-th down- 
and up-type quarks are proportional to the above Jarlskog invariants, which are of the 
ninth order in the Yukawa coupling, the quark (C)EDMs are highly suppressed in the SM. 

In the MSSM, the EDMs are not necessarily suppressed at the same level as in the 
SM. As long as one remains within the class of models belonging to the so-called minimal 
flavor violation (MFV) principle [22-24], where the flavor- violating terms in the MSSM 
arise only from the CKM matrix, the resultant EDMs are highly suppressed, at the same 
level of the SM. In contrast, within an MSSM framework with general flavor violation, 
the invariants introduced by the new flavor- violating interactions are less suppressed by 
the Yukawa couplings and large EDMs can be generated. 

Let us first parameterize the effects of flavor violation as usual by means of the mass 
insertion (MI) parameters [25,26], 

/ ™ \ _ {'^u)ij 
\"RR)ij — 2 ' 

(6) 

for i ^ j and rrij (/ = q, u, d, I, e) are the corresponding average sfermion masses. Notice 
that the above Mis are basis-dependent quantities. In this paper, we define the Mis in 
the bare SCKM basis, where the fermions and sfermions are rotated in the same way so 
that the corresponding tree-level Yukawa couplings are fiavor- diagonal; in contrast, the 
threshold corrections induced by non-holomorphic Yukawa couplings, will be off-diagonal 
in the ffavor space. We also take S'l^ = and the MI parameters for left-handed up- 

and down-type quarks are related by the SU{2)i symmetry as 5^j^ = VdiiV^ where V IS 
the hare CKM matrix appearing in tree-level Yukawa couplings. 

The size and the pattern of the Mis are unknown, unless we assume specific models. 
They are regulated by the SUSY-breaking mechanism and by the interactions of the high- 
energy theories beyond the MSSM. In this way, it makes sense to study the individual 
impact of different kinds of Mis on the low-energy observables. 

When only 7^ 0, the following JIs appear, 

jS) = Im{y,[K,5^Jh,, 

Jti^ - "^^{vAYa.SU},,. (7) 



(m 
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When {5'^j^)ij 7^ or 7^ 0, the new JIs are 

jJr^) = Im{5]^«y„F4,,. (8) 

Notice that both and jjj^^^ are of the third order in the Yukawa couphng constants 
and they arise from the relative phases between the corresponding Mis and CKM mixings. 

When both left- and right-handed squarks have mixings, the new JIs are proportional 
to only one Yukawa coupling constant as 

j'L:^ = l^{SRRyJL}u, (9) 

and the relative phases between the left- and right-handed squark mixings contribute to 
them. Notice that if only the left-handed squarks have mixing, the JIs are proportional 
to the external light quark masses. However, when the right-handed squarks have mixing, 
the JIs, and then the quark (C)EDMs, are enhanced by the Yukawa coupling constants 
of heaviest quarks. 

Finally, for the lepton EDMs, only the following JI is relevant, 

Jii = lm{6^j,j,y,6[^}^^. (10) 

This is because the leptonic sector does not have a counterpart of the CKM mixing within 
the MSSM. 

Here, we summarize the LO contributions to the flavored EDMs leaving the derivation 
of the BLO ones to the next section. The SUSY contributions to the quark (C)EDMs 
can arise at the one- loop level, and they are generated by the JIs of Eqs. ([7]) and ([9]) . In 
particular, when the left-handed squarks have mixing, the Higgsino diagram of Fig. ([I]) 
(on the right) contributes to the down-type quark (C)EDM through J^l^ and the resultant 
(C)EDM is given by 

(-• '''".] = ■;^^4^*. I- msiLhi {/f C(y)}, (11) 

where 2 = 1, 2, t/j = tan (3,y = /x^/m^ (m is the average squark mass) and {/5^'*(1), /o^''(l)} = 
{2/15, 1/10}, with the complete expressions for {f'^\y), fo'\y)} given in Appendix [El 
Eq. ffTT]) clearly shows that the (C)EDMs are suppressed by the external light quark 
masses, as expected from the invariants of Eq. ([7]). Assuming CKM-like Mis, i.e. isl ~ 

(0.2)^, a typical SUSY mass scale msusY = 500 GeV and tan/? = 10, d^/e and d'^ are of 
order ~ 10"^*^ cm. For up-type quarks, there is also an analogous expression for d^} 
proportional to J^^; however, it turns out that d^^^ /du} ~ [m^/(m^t^)] x [Ai/i//i^] x tp, 

(c) 

thus, dui is always negligible. 
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Figure 1: Leading contributions to the down-type quark (C)EDMs proportional to 
(left) and (right). 



At one-loop level, the are also gluino contributions to the quark (C)EDMs through 



the invariants J^j:^ (see Fig. ([T])). The induced down-type quark (C)EDMs are given as 



e 



t, Im [{Si.USUsr] {/f (^), /f (a^)} , (12) 



where i = 1,2, x = Ml/m"^ and {/f^(l), /f^(l)} = {4/135, 11/180}. Eq. ([HD shows 
that ddje and d"^. are enhanced by the Yukawa coupling of the heaviest quark and they 
can reach the level of ~ IQ-^^s-ae) cm for mgusY = 500 GeV, |(5^i)i3| ~ |(4/j)3i| ~ (0.2)^ 
and tan /3 = 10. Thus, the current experimental bounds on the hadronic EDMs already set 
constraints on the SUSY parameter space. Similarly, also the up-type quark (C)EDM is 
generated at the LO through a gluino mediated diagram and the corresponding prediction 
is 

where x = Mi/m'^ and {/f ^(1), /f ^(1)} = {-8/135, 11/180}. 

When only the right-handed squarks have mixing, the (C)EDMs are proportional to 
J^^l and they are generated only at the two-loop level. The most relevant effects are 
provided by non-holomorphic Yukawa coupling corrections, and they can still lead to 
sizable contributions, because i) they are enhanced by the heaviest Yukawa coupling via 
the CKM and the right-handed squark mixings, and ii) the non-holomorphic Yukawa 
couplings entering the (C)EDMs provide an additional tan/5 factor that can partially 
compensate the loop suppression e ~ as/(47r) if etan/3 ~ 1. 

Passing to the leptonic sector, the dominant contribution to the lepton EDMs arises 
from the one-loop exchange of binos/sleptons, and the corresponding EDM, proportional 
to J^^ , is given as 

4, ayrrirfiMi . , x /re x l i.(3)/ A r-,A\ 

~i = ^^^^^^ {{^ll\A^rr)^^ fo , (14) 
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where i = 1,2, x = Mf/m^ {m is the average slepton mass) and /q (1) = —1/15. The 
electron EDM approaches the level of ~ 10~^^ecm for msusY = 200 GeV, |((5i2.)i3l ~ 
Wrr)\ ~ (0-2)^ and tan/3 = 10. 

As stressed above, one peculiar feature of the flavored EDMs is that they can result 
to be proportional to the heaviest Yukawa couplings. This huge enhancement factor can 
bring the (C)EDMs close to the current and future experimental sensitivities, providing 
a splendid opportunity to probe the flavor structure of the MSSM. 

On the other hand, the predictions for the flavored EDMs are not enough accurate 
yet. In fact, in order to take into account the contributions arising from two- loop 

calculations are unavoidable. Moreover, for large values of tan j3, the LO predictions for 
the (C)EDMs discussed in this section receive significant corrections that have not been 
accounted for so far. In the following sections, we discuss the flavored EDMs at the BLO, 
stressing the importance of these contributions to get a correct prediction for the EDMs. 

3 Flavored EDMs at the BLO 

The effective Lagrangian necessary to evaluate all the relevant BLO effects to the (C)EDMs 
includes effective Higgs couplings with fermions and effective fermion-sfermion couplings 
with charginos and gluinos. The dominant BLO contributions are computed by including 
all the one-loop induced (tan /3-enhanced) non-holomorphic corrections for the charged 
Higgs (i?^) couplings with fermions and for the /g couplings with fermions-sfermions. 
The above effective couplings lead to the generation of if^ effects to the (C)EDMs, absent 
at the LO, via the one-loop i/^/top-quark exchange. Moreover, the chargino contribu- 
tions, suppressed at the LO by the light quark masses, are strongly enhanced at the BLO 
by the heaviest-quark Yukawa couplings. Finally, also the gluino effects receive large 
BLO contributions that are comparable, in many cases, to the LO ones. As a result, 
BLO effects dominate over the LO effects in large regions of the SUSY parameters space 
and their inclusion in the evaluation of the flavored EDMs is mandatory. The detailed 
implementation of the procedure leading to the above effective Lagrangian and the result 
in the mass insertion approximation are reported in Appendix [Al 

As stated before, the BLO expressions for the EDMs are obtained by inserting effec- 
tive (one-loop induced) vertices into the one-loop expressions for the EDMs. We would 
like to note that such an approach accounts for all the non-decoupling (tan /3-enhanced) 
contributions to the EDMs but it cannot provide the full set of two-loop effects. The 
latter requires a full diagrammatic calculation, which is outside the scope of this work. 

For instance, the expressions we find for the charged Higgs contributions will be 
valid as long as the typical supersymmetric scale msusY is sufficiently larger than the 
electroweak scale m^eak (~ t^Wi f^t) and the mass of the charged Higgs boson M^±. 
Therefore, our results can be regarded as the zeroth-order expansion in the parameters 
(m^g^]^, M^±)/mgusY of the full computation. However, it has been shown in Ref. [27] 
that this zeroth-order approximation works very well, at least in the 6 — >■ 57 case, even 
for Mh± > msusY, provided msusy is sufficiently heavier than m^eak- This finding should 
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also apply to our case, given that both b ^ s'y and the EDMs arise from a similar dipole 
transition. 

In addition, there exist potentially relevant two-loop effects, proportional to large log- 
arithms of the ratio wT.susY/''^weak, stemming from: i) the renormalizations of Yukawa cou- 
plings in the Higgs/Higgsino vertices, and ii) the anomalous dimensions of the magnetic 
and chromo-magnetic effective operators. These two classes of terms become important 
when the scale of the supersymmetric colored particles is significantly higher than the W 
boson and top quark masses. We account for these effects in the numerical analysis of 
Section HI 

However, we emphasize that the new two-loop effects we are dealing with in the present 
work are comparable to and often larger than the leading one-loop effects to the EDMs 
induced by gluino/squarks loops. Thus, a full inclusion of all the two-loop effects to the 
EDMs, although desirable, is not compulsory. 

Keeping the above considerations in mind, in the following sections, we derive the 
BLO formulae for the (C)EDM at the SUSY scale. The values attained by the (C)EDMs 
at the hadron scale (~ 1 GeV) are obtained according to Ref . [28] . 

For simplicity, we assume degenerate squark masses and we work in the mass insertion 
approximation, to make the dependences of the (C)EDMs on the SUSY parameters more 
transparent. 

Since the BLO effects to the quark (C)EDMs mostly affect the predictions for the 
down-type quark (C)EDMs, hereafter we do not discuss the subdominant BLO effects to 
the up-type quark (C)EDMs. 

3.1 Gluino Contribution to Down-type Quark (C)EDMs 

The relevant gluino-mediated contributions entering the evaluation of the down-type 
quark (C)EDMs at the BLO are reported in Fig. [2l As we can see, the leading BLO 
effects are basically obtained from the LO ones by means of the replacement of the tree- 
level couplings with the effective couplings derived in the Appendices. 

At the BLO, contributions proportional to all the JIs, J^^ , J^l' and J^^ naturally 
arise, in contrast to the LO case where only and J^l^ contribute to the (C)EDMs. As 
a result, the total BLO gluino-mediated contributions to the down-type quark (C)EDMs 
are given by 





(15) 
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dm djii bi dii d^i bi dn 




Figure 2: Gluino-mediated contributions for down-type quark (C)EDMs beyond the lead- 
ing order. In each diagrams gray boxes indicate the effective couphngs generated by 
non-holomorphic Yukawa corrections. 



where 



and 



' g 



1 + e3tp 3(1 + e^tp) V 1 + 63^/3 1 + 63^/3 

Wl + est.ni + est,) 6(1 + e,t,)^ J ^"^^ ' ^ 

(l + 63t,)(l + e3t,)^^ ^"^^ 3(l + e3t,)2(l + 63t,)^^ ^''^ 



' 4(l + e3t^)2 3(l + e3t/3)2(l + e3V 



+ 7 7/f(^), (18) 

(l + e3t/3)(l + e3t/3)^^ ^ ^' ^ ^ 



/f(^) = {-^/f(x), -lfi'\^) + lfi'\x)^ (^ = 1,2,3), (19) 

where x = Mf/m^ and the loop functions satisfy '^'^^(1) = {2/27, 5/18}, {-2/45, - 
7/60}, {4/135, 11/180} whereas their full expressions are listed in Appendix[El We also 
use the short hand notations e^, cr, clr, e^y, zy and zl instead of e^gj, eR^i, eiRu, eLyi, 
Re[zyJ and Re[zL.], respectively. The above loop factors are defined in Appendix [Bl and 
for equal SUSY masses and small MI parameters, they are given as eLR = €l = ^r = = 
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ei = sign(/iM3) x Os/Stt, ey = e^y = — sign(/iAt) x yf/32n'^, We also define 63 = £3 + ey 
and the rj and fj parameters such that 



1 + 63^/3 _ |l + eit;3| 



(20) 



Terms proportional to Sf, and contribute to the Jarlskog invariants , J^j^ and 
respectively. 

3.2 Chargino Contribution to Down-type Quark (C)EDMs 

As discussed in previous sections, chargino contributions to the down-type quark (C)EDMs 
arise already at the LO through a contribution proportional to J^i , hence suppressed by 
the light quark masses. At the BLO, new contributions proportional to both and 
J^^ (thus proportional to the heaviest down-quark mass mj,) arise. 

Combining the BLO contributions from the upper diagrams of Fig. [3] with the LO 
contribution of Eq. ffTTl) . the total charged Higgsino contribution reads 

+ £t^lm[ViML)3.]], (21) 

771}, I 



where 



- f (v) (22) 



and 



pi/± _ r-i (2) njeLY + eL)tp m) , . . . 



fgiy) = {lf!>'\y) - f?{yl f^\y) \ = i, 2) , (24) 



where y = ^"^/m^ and /]^±^(1) = {—5/18,-1/6}, {2/15, 1/10}. The charged Higgsino 
contributions play a very important role in the calculation of the EDMs resulting often 
dominant over the LO gluino-mediated contribution. Additionally, there exist also con- 
tributions from charged Higgsino-wino mixing diagrams as shown in the lower diagrams 
of Fig. [3l The final expression of all these contributions reads 



+ ^"Im [^sK^Ds.] + St^^'^lra m^iM I , (25) 

I fill 
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7, 9 




1, 9 



1, 9 




7- 9 



1, 9 



7> 9 




ti d-Li 





Figure 3: Chargino-mediated contributions for down-type quark (C)EDMs beyond the 
leading order. The effective couphngs with tan /3-enhanced corrections are depicted by 
gray boxes. 



where 



and 



Si 



^Rtl3 (1) 
^9) 



3(l + e3t/3)2^^±^± 



+ 



9(l + e3t^)2(l + e3t/3) 6(l + e3t;3 



2 I ^H±W± 



[1 + ri)eReYt 



3(1+ 63^/3)2(1 + e3t/3)''^*^ 



/3 



3(1 + e^t^y 3(1 + -e^tpY{l + e^tp) J 



,(0) 



+ 



,(1) 



[x,y) 



9}, 



^^^w^^^^y") = \^3o\x,y) - 9?{x,y) , gl'\x,y)\ (i = 0, 1) 



(26) 
(27) 



(28) 



(29) 



where x = Ml/m^, y = fi^/m^ and gfl]-^{l,l) = {-11/18, - 1/6}, {13/45, 2/15}. 



The exphcit expressions of the loop functions are presented in Appendix [El 



3.3 Neutralino Contribution to Down-type Quark (C)EDMs 

Neutralino-mediated contributions are divided into three classes: pure bino-mediated dia- 
grams, neutral wino-Higgsino mixing diagrams and bino-Higgsino diagrams. The relevant 
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diagrams for the neutralino contributions are reported in Fig. |H The pure bino-mediated 
contributions are given by 

e J ^ 47r \ '- -' 

+ El Im + El ^Im {V^S^^iM ) ' (30) 



where 



and 



/ (1 + n)e^entl _ r^e^Btp \ (i) . . 

^ 1^9(1 + e3t;3)2(l + e3t/3) 6(l + 63t;3)V^ ^ ^' ^ ^ 

^" (l + 63t,)(l + 63t,)^^ ^"^^ 3(l+63t,)2(l + 63t,)^^ ^^'^ 

^ V3(l + e3t;3P 3(l + 63t^)2(l + e3t/5)r^ ^ ^ 



/^'^(^) = <! -\fo\^), (^) ^ = 1,2,3) , (34) 



3' 



where x = M^/m^ and 4^'^'^Hl) = {1/18, -1/6}, {-1/30, 1/10}, {1/45, -1/15}, respec- 
tively. Here, the hypercharges and Yr are defined as = 1/6 and Yr = 1/3. 
Next, the contributions of neutral wino-Higgsino mixing diagrams read 

+ £:f ^°Im [V*Mn)s.] + El'^'^lra [V*{5i,u] , (35) 



where 



3(1 + £3^/3) 



9(1 + e3t^)2(l + e,tp) 6(1 + e3t;3) V 



f2Jwo(a;'2/) ' (36) 
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.(1) 



9'^Loi^,y), (38) 



and 



9%^oi^,y) = [-l9iH^,y) , 9^\^,y) \ = o, i) , (39) 

where x = Ml/m\ y = i^^m^ and gf^l^il, 1) = {1/18, -1/6}, {-2/45, 2/15}. 
Finally, the contributions from bino-Higgsino mixing diagrams read 



HB 



+ ^Im [K,l(4ii)3.] + ^a^^^Ini [V^3:(<^^l)3.] ) , (40) 



where 



nHB_ I 'h'^.i y , y \ (1) / N 

+ f - .7?^) + y) , (41) 



9(1 + 63^/3)^(1 + 63^/?) 6(1 + 63%; 



sAb_ Ji) . . _ (yL + yi.)(l+r.)e^eyt^ 

- (l + 63t,)(l + 63^,)^^^^'^'^^ S{l + est,ni + e,t,) ^^B^'''y>' ^^'^ 



3 V3(l + e3t/3)' 3(l + e3t^)2(l + e3i^); 



(0) / N 

.3(1 + e.tpf 3(1 + e3t^)2(l + 63^^) J 

^ (l + 63t,)(l+63t,)^^^^^'^^' ^^'^ 



and 



4y^'^) = ^-^5^^(^'^)' ^^^(^'^)^ (^ = 0,1), (44) 

where x = Ml/m?,y = iJ?/m? and ^|^^(1,1) = {1/18,-1/6}, {-2/45,2/15}, respec- 
tively. 
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^Ri bn bi dju dm d^ d^ dj^ bi d^ dn 




Figure 4: Neutralino-mediated contributions for down-type quark (C)EDMs beyond the 
leading order. Gray boxes indicate the effective couphngs with possible loop factors. 
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Figure 5: Left (Right): Charged Higgs mediated contributions to down-type quark 
(C)EDMs at the leading order for the case of J^^ (Jlr)- The effective couphngs with 
tan /3-enhanced corrections are depicted by gray boxes. 



3.4 Higgs Contribution to Down-type Quark (C)EDMs 

The importance of the charged-Higgs mediated contributions to the hadronic (C)EDMs 
was aheady discussed in Ref. [12] where, on the other hand, only the leading contribu- 
tions of order Oletp) were presented. The full expression for the charged Higgs-mediated 
contributions, beyond the leading-order expansion in Oi^etp), reads 



e 



a2 mb 

16^ mJI^ 



+ St Im [v;MR)3^] + —£t Im [viMM 

mi) 



(45) 



where 



and 



^Rtp 



9(1 + est^) 

^Rtp 



1 + 63^/3 



3(1 + 63^/3)2 



1 + 63^/3 



(0) , 



+ 



3(1 + e3t/3)(l + 63^/3) 3(1 + 63^/3 



1-^3^/3 + 



'?>{l + e;tp){l + e^tp) 



4^/3 



/3 



1 + 63^/3 



1 + ^3^/3 



/ 



(0) , 



(46) 
(47) 



(48) 



/o 



(0), 



(49) 



where z = mf/Mfj+ and /^i(l) 



{—7/9,-2/3}. The loop factors with primes come 
from the vertex correction for charged Higgs with left-handed down quarks, and these 
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are given as = £3 = siga^fiMs) x as/Svr, e'y = —sign{fiAh) x ?/^/327r^ for equal SUSY 
masses. 

We note that, since the loop function f['^\z) contains a large logarithmic and it behaves 
as f['^\z) ~ 3 + logz for z -C 1, both EDMs and CEDMs have a slow decoupling with 
the charged Higgs masses going as d^^^ ~ z log z. Even for moderate tan f3 values and 
charged Higgs masses Mfj± < 2 TeV, the down quark (C)EDMs remain within the reach of 
sensitivities of future planned experiments, i.e. {dd/e)H± and {d^)H± ~ C(10~'^^^~^^)) cm 
for I (5^/;;) 13 1 = (0.2)^. Another important feature of the charged Higgs contribution is 
that it does not decouple for msusY — > 00, as long as the Higgs mass doesn't decouple. 

Finally, additional contributions to the (C)EDMs come from the neutral-Higgs sector. 
However, the H and A contributions have opposite signs and their net effect to the 
(C)EDMs is really small even in the most favorable case where a mass splitting between 
H and A of order 10% was allowed. 

Before concluding this section, we would like to note that the flavored (C)EDMs are 
highly correlated with FCNC observables. As an interesting example, let us mention that, 
irrespective to the particular choice for the SUSY spectrum, the and contributions 
to the EDMs are related to the NP contributions entering B Xs'y as 

+ {d,)H. ^ -eg^Txrr [^3:(4ij)3.] c, , (so) 

where Cj = C^^{Mw) + C^^{Mw) is defined as B{B^X,j) ~ 3.15 - 8 C7 - 1.9 Cg [29]. 
A detailed exploration of the intriguing correlation and interplay among FCNC processes 
and flavored (C)EDMs would deserve a dedicated study that goes beyond the scope of 
the present work. 



3.5 Comparison of Various Contributions for Down-type Quark 
(C)EDMs 

As we have discussed in the previous sections, the down-type quark (C)EDMs receive the 
first contributions already at the LO level through the exchange of gluinos and charginos. 
Moreover, at the BLO, many additional contributions, proportional to the various JIs, are 
generated. The purpose of this section is to provide a detailed numerical comparison of 
the various effects. In the following, we discuss the contributions from and J^j:^ which 
are typically more important than the contributions from J^l^ thanks to the enhancement 
by the heaviest Yukawa couplings. 

In Fig. [6] (a) [(c)] and (b) [(d)], we show various contributions to the down quark EDM 
[CEDM] as generated by the phases 0i = Im [(5^^)13(4^)31] and 02 = Im [■i^3*i(4i?)3i] • 
We choose CKM-like mixing angles |(4l)3iI — l(4i?)3il = (0.2)'^ and we assume maxi- 
mum phases; moreover, we consider a common soft SUSY mass msusY = 1 TeV, Mh± = 
500 GeV and a positive fi term. 

As shown in Fig. [6] (a) and (c), the gluino effects are the dominant ones as long as 0i 
is the only non zero phase. We note that the LO gluino effects have a linear dependence 
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dd (e cm 




Figure 6: Various contributions to the down quark EDM in (a), (b) and CEDM in (c), 
(d) as functions of tan j3. Here, the masses of SUSY particles arc taken to be ITeV while 
Mh± — 500 GeV. Contributions from the phase 0i = arg[(5|^^)i3((5^^)3i] is presented in 
(a) and (c) and those from the phase 02 = ^'^zWiii^Rid^^ ™ (b) ^^<i (d)- ■'■^ these 
figures, maximum phases are assumed and only the absolute values of each contribution 
are plotted. 
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with tan /?, however, the inclusion of higher order corrections acts to reduce sizably the 
LO effects, specially in the large tan (3 regime. 

The wino-Higgsino mixing contributions turn out to provide the most important effects 
after the gluino ones. In particular, they almost grow as tan^ (3 and become comparable to 
the gluino effects at around tan (3 ~ 60. However, in realistic mass spectra, the charginos 
tend to be lighter than the gluino, thus, the wino-Higgsino contribution becomes more 
important compared to what Fig. [6] (a) and Fig. [6] (c) show. 

Turning to the EDMs from 02, which are generated only at the BLO, we find that the 
contributions from charged Higgs and Higgsino are dominant for any value of tan/5. In 
particular, for low to moderate values of tan /3, the charged Higgs contributions (growing 
almost linearly with tan (3) provide the dominant effects while from moderate to large val- 
ues of tan/3, the Higgsino contributions dominate thanks to their quadratic dependence 
on tan [3. Fig. [6] (b) and Fig. [6] (d) also show that the linear dependence on tan (3 of the 
charged Higgs effects is partially lost at large tan /5 due to higher order effects. Contribu- 
tions from the gluino are typically subdominant but still not completely negligible. 

Next, we compare charged Higgs and Higgsino contributions from 02 with the gluino 
contribution from 0i. Their relative size is controlled by the ratio | V'3\/(5^^)i3|, which we 
set to be 1 in Fig. [61 In concrete models, as for instance the CMSSM or SUSY 5*^7(5) 
models with right handed neutrinos, is generated by radiative correction through the 
CKM matrix and the top Yukawa coupling, and it turns out typically that (5^2,) 13 — —cV^^ 
with c ~ (9(0.1). As a result, the BLO charged Higgs and Higgsino contributions can be 
dominant over the LO gluino contributions. 

All the above considerations for the EDMs apply also to the CEDMs with the only 
difference being that the chargino contributions for the CEDMs are less important than the 
EDM ones because of the absence of some constructive interferences among amplitudes. 

3.6 CP-violating four-Fermi interactions 

It is well-known that CP- violating four- Fermi interactions can generate atomic and hadronic 
EDMs through the effective (CP-violating) couplings of the neutral Higgs bosons with 
fermions [17,30,31]. These contributions to the EDMs are proportional to tan'^/3 and 
become important only for large tan (3 values and for Higgs masses much lighter than the 
soft masses. In Refs. [17,30,31], flavor-conserving but CP- violating phases were assumed. 

As we will see, CP-violating four-Fermi interactions can be also generated by fla- 
vor effects. In fact, as discussed in the Appendices, flavor-mixing effects contribute to 
both flavor- violating and flavor-conserving CP- violating couplings of the Higgses with the 
fermions (see Eq. (11061) ). 

The relevant CP-odd four-Fermi interactions generated by flavor effects have the fol- 
lowing form 





'3 ' 



(51) 



ft,fj=e,d,s,b 



with the coefficients C 



fifj 



deflned as 



19 



Im 



m , 



■h 



33 

Af--A 



fi 



4mlr 



m. 



\l + eitp + efhl\\l + e,tp + efhV 



(52) 



^3 "Z?! 



^^i*^^ given in Appendix [Xj is the effective coupling of the neutral Higgs H 
with the fermion / and Af. is such that 



A 



nib 



fi 



[1 + ri)eLeRtl 



L6(l + e3t^) 9(1 + 63^/3) (l + est^ 



[1 + ri)eReYtl 



rrid^ 3(1 + e3t^)(l + e-^tf, 



-Im [V;,{5i^),,] 



(53) 



for fi = d,s and ~ 0. If = e, the following expression arises 



A. = ^ 



rUe 



L6(l + e^t^) 9(l + e^t;3)(l + e^t 



/3JJ 



Im [((5; 



LL 



(54) 



where the loop factors e^, and e^^ are defined in Appendix 



masses it turns out that e. 



and 



= el = -sign(/iM2) x 3^2/167?, e^^ 



and for equal SUSY 
-- sign (/i Ml) X ay /Sit 



0; moreover, - ^ ^3.^;/ |x ^ c^.^ ^ .^l 
The most prominent effect of the above CP-odd four-Fermi interaction compared to 
the flavor conserving case of Refs. [17,30,31] is that, thanks to the flavor effects, we can 
always pick up the heaviest Yukawa couplings in the Higgs couplings with fermions. Let 
us note that in Eq. fl53l) we have neglected the contributions proportional to J'll because 
they are not enhanced by the heaviest Yukawa coupling, hence, they can never reach an 
experimentally interesting value. 

In order to understand the impact of the flavored four- Fermi interactions on the atomic 
and hadronic EDMs, it is convenient to consider the two Higgs doublet model (2HDM) 
limit of the MSSM, where the soft sector is assumed to be much heavier than the Higgs 
sector. In this limit, de is negligible and (ixi can be only generated by the CP-odd four- 
Fermi interactions (see Eq. ([2])) while the main contribution to dn arises from the charged- 
Higgs mediated diagrams. In Fig. [TJ we show the predictions for d^i and dn in the 2HDM 
limit of the MSSM setting Ag = and assuming a common heavy SUSY mass for the 
soft sector. Here, we assume that the EDMs are generated by the flavor dependent 
CP phases in the squark mass matrices. The shaded region in Fig. [7] is excluded by 
the current experimental bound on dn- Considering the current experimental bounds on 
d^i < 9 X 10~^^e cm [7] and dn < 3 x 10"^^e cm [6], we conclude that dn represents 
the most powerful probe of the 2HDM limit scenario, when CP violation is generated by 
flavor effects in squark mass matrices. 
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Figure 7: (ixi vs (i„ in the 2HDM limit of the MSSM. Here, we assume that the EDMs are 
generated by the flavor dependent CP phases in the squark mass matrices. The shaded 
region is excluded by the experimental bound on dn 



3.7 Charged lepton EDMs 

As already discussed in Section [2|, the dominant contributions to the lepton EDMs arise 
at one- loop level through the invariant J^^ , that is generated by the exchange of binos 
and sleptons. 

On general ground, BLO corrections to the charged lepton EDMs are not so important 
as for the down quarks C3iSG, clS the gluino and Higgsino-mediated contributions are absent 
in the leptonic sector. 

In the following, we present the complete formulae for the charged lepton EDMs at 
the BLO, which include the exchange of bino/sleptons, bino-Higgsino/sleptons and wino- 
Higgsino / sneutrinos 



9 9 t/J 



-HB 



ay Ml 



Im[(5i,).3(5|«)3.] , (55) 



where 



, (4 + ^ytp A2) 



(1), , 
6(i + e^t^)2^o ^^i.y) + 



Xl) + 



+ r^)ele'ktl 



[12(1 + eltpy 18il + eltp) 



^1 



HW 



R'-P 



9(1 + emf 



9frJx2,y), (56) 



and the loop functions for wino-Higgsino mixing diagrams are given as 



9^f!w(^2,y) = -gi'{x2,y) - ^g^''{x2,y) {i = 0,1) 



(57) 
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Figure 8: d^^*^ /d]f^ vs tan/? taking the /i term, the slepton and the gaugino masses 
randomly in the range of [200 : 1000] GeV. Light gray points satisfy the GUT relation for 
the gaugino masses, M2 = 2Mi, while dark gray points do not. 



where X2 = M|/m^, Xi = M^/rn^ and y = fi'^/m'^, with m being a common slep- 
ton/sneutrino mass. 

In Fig. [HI we plot the leptonic EDMs evaluated at the BLO normalized to the LO 
estimate as a function of tan f3, scanning over the relevant SUSY masses in the range [200 : 
1000] GeV. The inner region of Fig. [8] (light gray) corresponds to the points satisfying 
the GUT relation between the gaugino masses, M2 = 2Mi, while the points in the larger 
region (dark gray) do not satisfy such a GUT relation. 

As shown by Fig. [HI BLO corrections to the leptonic EDMs can provide only modest 
contributions, at the level of (9(10)% for large values of tan/3, if the GUT relation between 
the gaugino masses is assumed. In contrast, if we do not impose such a GUT relation, large 
BLO effects can arise in the large tan/? regime at the level of 0.5 < d^^'^/d]:^ < 1.5. Such 
a behavior of d^^'^/d];^ with the gaugino masses can be understood noting the different 
dependence of d^^'~' and d^.^ on the gaugino masses; in particular, when M2 <^ Mi, it 
turns out that the LO contributions have a faster decoupling property with Mi compared 
to the BLO contributions. 



4 Flavored EDMs in SUSY GUT models 

In this section, we discuss the predictions for the (C)EDMs in some specific SUSY models 
as SUSY SU{5) with right-handed neutrinos. As we will show, in this class of models, 
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the inclusion of the BLO contributions presented in this work turn out to be crucial to 
get the correct predictions for the (C)EDMs. 

Since SUSY GUT models present rich flavor structures, they typically predict non- 
negligible flavored (C)EDMs [32] as well as non-negligible rates for various flavor- violating 
phenomena [33]. Moreover, since within SUSY GUTs leptons and quarks sit into same 
multiplets, the flavor violation in the squark and slepton sectors may be correlated [33], 
if the SUSY-breaking effects are mediated at such a high energy scale at least. 

A detailed exploration of the intriguing correlation and interplay among FCNC pro- 
cesses in the leptonic and hadronic sectors with the flavored (C)EDMs would be very 
desirable, but this analysis would require a dedicated study that goes beyond the scope 
of the present work. On the contrary, the major aim of this section is to stress the impor- 
tance of BLO effects for the flavored EDMs in well motivated SUSY scenarios as SUSY 
SU{5) with right-handed neutrinos. 

Let us first review the flavor structures in the squark and slepton mass matrices in 
the SUSY 5*^7(5) GUT with right-handed neutrinos. The Yukawa interactions for quarks 
and leptons and the Majorana mass terms for the right-handed neutrinos are given by the 
following superpotential, 

W = ^yu,,^i^jH + V2yd,j^,^,Tl + y,i^l^i^,H + Mi,l^,l^,, (58) 

where \E' and $ stand for 10- and 5-dimensional multiplets, respectively, and refers 
to the right-handed neutrinos. H {H) is a 5-(5-) dimensional Higgs multiplet. After 
removing the unphysical degrees of freedom, the Yukawa couplings and the right-handed 
neutrino mass matrix of Eq. fl58l) can be written as 

\ / ij 

Vdij = VdAj^ 

yu,, = e-'^-^W^e-'^-^yM\ 

Mij = M,Aj. (59) 

In this paper, hat and bar symbols refer to diagonal matrices and to bare quantities, 
respectively. In Eq. (I59l) . (pd, (pu and 0p are the physical CP- violating phases satisfying 
(j) f. = {f = u, d, V and v). The unitary matrix V is the CKM matrix in the extension 
of the SM to the SUSY SU (5) GUT, and each unitary matrices \J, V and W have only a 
phase. In the limit where W = 1, which we assume, U corresponds to the MNS matrix. 
In this case, the light neutrino mass eigenvalues are given as 

m.. = ^(Hj)\ (60) 

where Hf is a doublet Higgs in H. The colored Higgs multiplets He and He are introduced 
in H and H respectively as SU (5) partners of the Higgs doublets of the MSSM. They have 
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new flavor- violating interactions in Eq. (158!) . If the SUSY-breaking terms in the MSSM are 
generated by dynamics above the colored Higgs masses, such as in the gravity mediation, 
the squark mass terms may get sizable corrections by the colored Higgs interactions. In 
the mSUGRA scenario, low-energy flavor-violating SUSY-breaking terms are radiatively 
induced, and they are qualitatively given as 



[m 



2 X 2ta\ /o 2 , a2\ / o t , 1 -^GUT \ 



iVmVU^ml + Al) (3 In -Fin. 

V Mgut 



i^ny' ;^n-'-o-^^o; Mgut ^ msusY; 



= -(ip(^'^^«^*)^^(3-o + ^g)ln^, (61) 

with 2 7^ j; we have also assumed the colored Higgs mass to be the GUT scale Mqut and 
mo and Aq are the universal scalar mass and trilinear couplings, respectively As shown 
by Eqs. ( pTl) . the off-diagonal components of the right-handed squark and left-handed 
slepton mass matrices are induced by the neutrino Yukawa couplings, and they depend 
on the various CP- violating phases of the model. The existence of flavor violation in the 
right-handed down-type squarks is essential to generate the contributions to the down- 
type quark (C)EDMs proportional to the Jarlskog invariants J^j^^ and J^^^, which are 
enhanced by the largest down quark mass m{,. In this stressed in the previous 

section, not only the gluino but also chargino and charged Higgs are expected to contribute 
significantly to the (C)EDMs. 

Although in the previous section we have reported the analytical expressions for the 
down-type quark (C)EDMs in the MI approximation, our numerical analysis is performed 
using the mass eigenstates. Moreover, we also include the leading QCD corrections and 
we derive the (C)EDMs at the hadron scale (~ 1 GeV) [28]. 

In Fig. [9l we show the various contributions to the down quark (C)EDM, setting mo = 
Ml/2 = 400 GeV, Aq = and /i > 0. For the neutrino sectors, we assume a hierarchical 
mass structure for light neutrino masses and take rrii,^ = O.OSeV, Mj^g = 10^^ GeV and 

Ue3 = 0.01. 

For the down quark EDM in Fig. [H] (left), the chargino and charged Higgs contributions 
are typically comparable to the gluino ones despite the fact that they are generated 
only at the BLO. This happens because i) the charged Higgs, chargino and stop masses 
(entering the BLO contributions) are generally lighter than the gluino and the other 
squark masses (entering the LO contributions) in most of the MSSM parameter space; ii) 
\{^ll)i3\ < 1^1 1 5 when is radiatively- induced via the CKM matrix (as it happens in 
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Figure 9: The down quark EDM (left) and CEDM (right) in the SU{5) GUT model 
with right-handed neutrinos as functions of tan/5. Various contributions to the down 
quark EDM and CEDM are also shown. Here, maximum phases are assumed and only 
the absolute values of each contribution are plotted. The input parameters are given as 
mo = Ml/2 = 400 GeV, Aq = and /i > 0. The predicted values are approximately 
proportional to the mass of right-handed tau neutrino and the (1,3) component of the 
MNS matrix, which are taken to be M^,^ = 10^^ GeV and 11^3 = 0.01, respectively. The 
region for tan/? > 35 is already excluded by the current experimental limit on the LEV 
process, Br{^ — > e-y) < 1.2 x 10^^"*^. 



our model); in) the mass functions for the chargino (charged Higgs) contributions tend 
to become larger than those for the gluino contributions due to positive interference of 
diagrams (due to the existence of large logarithm). 

Moreover, when tan/? is large, the BLO effects become more significant and they 
dominate over the LO ones. In the chargino case, this is explained by the explicit tan j3- 
dependence of {dd)x^± ~ t^, to be compared with {dd)g ~ t/?. In the charged Higgs case, in 
spite of the same exphcit tan /3-dependence of {dd)g and {ddjn, {d^jH > {dd)g for increasing 
tan/?, since uih is highly reduced by large RG effects driven by ^ ^ 1. 

Notice that, the corner of the MSSM parameter space where the BLO effects 
are particularly enhanced compared to the LO ones, corresponds to the so-called A- 
funnel region (where tua — 2mLSp), satisfying the WMAP constraints. In the mSUGRA 
scenario, for yU > 0, the chargino contributions have an opposite sign compared to the 
gluino and charged Higgs contributions, leading to a large cancellation between them. 

From Fig. M we note that the chargino contribution to the CEDM is much lesser 
than the corresponding contribution to the EDM; this is because the EDM has positive 
interferences among diagrams that are absent in the CEDM case. 

The predictions for the down quark EDM are investigated in more detail in Fig. [TOl 
In particular. Figs, (a), (b) and (c) describe the value of the EDM in the (M1/2, mo) plane 
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for /i > and tan/? = 10,25 and 50, respectively. Fig. (d) refers to the case in which 
tan/3 = 25 and /i < 0. 

For relatively small values of tan/3 (see Fig. [TOT a)). the chargino contribution is not 
very large compared to the gluino and charged Higgs ones. On the other hand, while in- 
creasing tan /5, the chargino contributions become more important and they dominate over 
the other contributions in many regions of the parameter space, see Fig. [TOlfb). Chargino 
and charged Higgs contributions are specially important in the region where M1/2 > ^o- 
Although the off-diagonal components of the squark mass matrices mainly depend on the 
values of mo and Aq, as can be seen in Eq. ( 16T1) . their diagonal components are mainly 
determined by M1/2 due to the large RG effects driven by the SU{3)c interactions. As a 
result, the MI parameters become small for M1/2 > "^0 and the gluino contribution pro- 
portional to Iiii[((5^£,)i3(5^^)3i] is more suppressed than the chargino and charged Higgs 
contributions that are proportional to Im[l3*]^ (5^^)31]. 

This feature is clearly shown in Fig. [TT] which describes the ratio of the chargino 
and charged Higgs contributions over the gluino ones. In the region where M1/2 S> rriQ, 
dd{x^)/dd{g) approximately approaches to a constant for a fixed value of uiq. This is 
because, as mo/Mi/2 becomes small enough to be neglected, all the SUSY-breaking terms, 
even the off-diagonal ones, are mostly determined by M1/2, and all the contributions to the 
EDM behave as dd oc l/Mlj^. On the contrary, the charged Higgs contribution decouples 
more slowly due to the large logarithm in the loop function, thus it provides the dominant 
effects in the regions with heavy SUSY particle. This parameter dependence explains the 
features of Fig. [TOT b): for small M1/2, the gluino is relatively light and it provides the 
most important contribution to the (C)EDMs while, as M1/2 increases, the charged Higgs 
and specially the chargino contributions start dominating. For very large M1/2 values, 
strong cancellations among amplitudes do not appear because of the dominance of the 
charged Higgs contribution. 

For very large tan/5 (~ 50), the BLO chargino mediated contributions tend to become 
the dominant one in large regions of the parameter space (see Fig. [TOlfd)) due to their 
tan^/3 dependence. In the case of /i < (see in Fig. [TOlfd)). the down quark (C)EDMs 
reach values significantly larger than those for /x > 0. In fact, when /i < 0, gluino, chargino 
and charged Higgs contributions have a common sign, thus, they constructively interfere 
with each other; moreover the threshold corrections to the tree-level Yukawa couplings 
provide enhancement contributions, in contrast to the /i > case. 

In summary, the inclusion of BLO effects for the (C)EDMs is mandatory as long as 
they are generated by flavor violating terms in the MSSM soft sector. In addition, when 
large cancellations among various contributions do not occur, the predicted values for the 
down quark EDMs can be within the reach of the future experimental sensitivities such 
as for the neutron EDM {dn ~ 10~^^ ecm) and deuteron EDM {dd ~ 10~^^ ecm). So, the 
observation of hadronic EDMs at these future experiments could allow us to explore the 
flavor- violating soft sector of the MSSM. In this respect, an accurate theoretical prediction 
for the EDMs, by means of the inclusion of the BLO effects discussed in this work, plays 
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Figure 10: Predictions of the down quark EDM in SUSY SU{h) GUT model with right- 
handed neutrinos in the (mo, M1/2) plane. The input parameters mo, M1/2 and Aq are 
given at the Planck scale and Aq is set to Aq = 0. The sign of the fi term is taken always 
positive but in (d) and tan/? is taken to be 10, 25, 50 and 25 for (a), (b), (c) and (d), 
respectively. The predicted values for the EDM are approximately proportional to the 
mass of right-handed tau neutrino and to the (1,3) component of the MNS matrix, which 
are taken to be Mjvg = 10^^ GeV and Ues — 0.01, respectively. The grey region is excluded 
by the current experimental bound Br{ii — > 67) < 1.2 x 10~^^. 
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Figure 11: Left(Right): Ratio of the charged Higgs-mediated (chargino-mediated) contri- 
butions over the gluino- mediated contributions in SUSY SU{5) GUT model with right- 
handed neutrinos. The input parameters are the same as in Fig. [TUT b). 



a crucial role B 

For completeness, we also discuss the predictions for the electron (muon) EDM de 
(d^), within the SUSY SU{5) model with right-handed neutrinos. In Fig. [121 we show 
the values attained by d^ and in the (mo, M1/2) plane setting tan/3 = 25 (for different 
values of tan/3, de and d^ approximately scale as (tan/5/25)). The grey region is excluded 
by the current experimental bound Br{^ 67) < 1.2 x 10^^^. As we can see, de and 
dfj_ can reach values at the level of de ^ 10~^^ ecm and df^ < 10^^^ ecm, well within the 
expected future experimental sensitivities on de, at least. 



5 Conclusions 

The hadronic and leptonic electric dipole moments (EDMs) represent an ideal ground 
where to look for New Physics effects. In fact, on the one side, the Standard Model 
predictions for the EDMs are well far from the present experimental resolutions and, on 
the other side, the EDMs are among the most sensitive observables to New Physics effects 
in a broad class of theories beyond the Standard Model, included supersymmetry. 

^ The hadronic EDMs could be also generated by flavor-conserving but CP-violating parameters, such 
as the /i term and the A terms. Actually, if non-vanishing EDMs will be observed in some experiments, 
powerful tools in shedding light on the source of CP violation would be correlated analyses of the 
EDMs for many kind of nuclei and atoms, and also ii) studies of the correlations among EDMs and flavor 
violating transitions [34] . Obviously, a precise calculation of the EDMs is essential for this purpose, too. 
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Figure 12: Predictions for the electron (muon) EDM (d^) in SUSY SU (5) with right- 
handed neutrinos in the (mo, M1/2) plane. The values attained by de and (i^ scale with 
tan/3 as (tan/5/25). See the caption of Fig. [TU]for additional details. 



However, there are two main obstacles to fully exploit the New Physics sensitivity of 
the EDMs: 

• Experimentally, one measures the EDMs of composite systems, as heavy atoms, 
molecules or the neutron EDM while the theoretical predictions are relative to the 
EDMs of constituent particles, i.e. quarks and leptons, thus a matching between 
quarks and leptons EDMs into physical EDMs is necessary and this induces sizable 
uncertainties related to QCD and nuclear and atomic interactions. 

• From a theoretical side, the supersymmetric evaluation of the quark and lepton 
EDMs has been performed only at the leading order, thus, potentially large effects 
from BLO contributions have been disregarded so far. 

In this paper, we have reported a detailed analysis, within supersymmetric theories, of 
the hadronic and leptonic flavored EDMs, i. e. EDMs generated by flavor- violating soft 
terms, at the BLO contributions, providing an answer to the second question. 

As shown by the present analysis, the impact of BLO corrections to the flavored 
EDMs, when arising from (tan /^-enhanced) non-holomorphic contributions, can be very 
important. We have pointed out the existence of new contributions to the EDMs that 
can arise only at the two-loop level, thanks to the inclusion of (tan /3-enhanced) non- 
holomorphic corrections. In fact, at one-loop level, certain contributions to the EDMs are 
absent (as it happens for the charged-Higgs mediated EDMs) or completely negligible (as 
for the chargino mediated EDMs). In the above cases, if new sources of flavor violation are 
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present, two-loop effects do not represent just sub- leading corrections to the corresponding 
one-loop effects but, on the contrary, they provide the largely dominant effects. 

Moreover, the most important result of our work is that these two-loop effects (from 
charged-Higgs and chargino exchanges, respectively), previously neglected in the liter- 
ature, are systematically comparable to and often larger than the one-loop effects from 
gluino exchanges in a broad region of the SUSY parameter space. This result is in contrast 
to what typically happens in flavor-changing neutral current transitions. For instance, in 
the 6 — ^ 57 transition, the leading contributions from charged-Higgs and charginos arise 
already at the one-loop level and BLO effects can only provide a sub-leading correction to 
the leading order result. Moreover, altough BLO effects to the EDMs arise from (tan/3- 
enhanced) non-holomorphic contributions, they turn out to be very important for the 
entire allowed range of tan /3. 

As a result, the evaluation of the BLO effects to the (C)EDMs presented in this work is 
an essential step forward towards a correct prediction for the (C)EDMs in supersymmetric 
theories with flavor-violating soft terms. This point has been explicitly demonstrated in 
Section m for the case of a SUSY SU{5) theory with right-handed neutrinos. 
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A Effective Interactions beyond the Leading Order 



The impact of non-holomorphic (tan /^-enhanced) Yukawa couplings has been aheady 
discussed in various processes [23,35-37]. However, in spite of their importance, a sys- 
tematic analysis including general flavor structures and CP-violating phases has not been 
accomplished yet. In this section we derive various effective couplings induced by the 
tan /^-enhanced radiative corrections with explicit CP phase dependence. 

The derivation of the effective couplings is performed by means of the following proce- 
dure [23,36]: i) evaluation of the one-loop induced coefficients for the dimension-four op- 
erators, which modify the Yukawa couplings, in the MI approximations, ii) determination 
of the unitary matrices that diagonalize perturbatively the quark and lepton mass matri- 
ces, treating the flavor-violating one-loop corrections as small parameters. At this step, 
we specify the relation between the tree-level (bare) couplings, that appear in the tree- 
level Lagrangian, and the observed (physical) parameters, in) derivation of the relevant 
effective couplings in the basis where the quarks and leptons are in the mass eigenstates 
while the squarks and sleptons are in the flavor eigenstates. 

In the next subsection we address the points i) and ii) while in the last two subsections 
we focus on the point Hi). 



A.l Non-holomorphic Yukawa CoupUngs 



In the SU{2)l x U{1)y symmetric limit, the interaction Lagrangian for the Higgs and 
matter fields is expressed by 



C 



if.. 



HLj 



+ 4. 



/°^. + h.c.. 



(62) 



where V and y^ are CKM matrix and Yukawa couplings, respectively, in the hare CKM 
basis and e{j (/ = u, (i, e) are the non-holomorphic radiative terms, which are gener- 
ated when heavy SUSY particles are integrated out from the effective theory. Here, the 
conjugate fields are defined as 

Hl^ ^ eatiH^.r, Hi ^ ea^iHuY • (63) 

The non-holomorphic radiative corrections are generated at one-loop level and they are 
reported in AppendixO Notice that Eq. (jHSD contains both the corrections for the fermion 
mass matrices and the vertex corrections for Higgs couplings. 

After the electroweak symmetry breaking, we obtain the fermion mass matrices which 
are diagonalized by the unitary matrices through 



"72 
'V2 



Vei {^^3 - tp^j) = iPl/rieU, 



(64) 
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where sp = sin/?, = cos/3, tp = sin/3/ cos/? and v = 2mw 1 92 while rhf. stands for the 
mass eigenvalues with Uf (with / = ul, ur, (11, dn, cl, cr) being the unitary matrices that 
diagonalize the quark and lepton mass matrices. In the above equations, hat and bar 
symbols refer to diagonal matrices and hare quantities, respectively. The physical SCKM 
basis is obtained by the unitary transformations 

Ul = Uu^V ul, dL = Ud^dl, cl = Ue^el, fn = Uf^^f^ (/ = u, d, e) , (65) 

and the CKM matrix is expressed by 

V = UulVUI . (66) 

In order to include the above threshold corrections to the Yukawa couplings, we need to 
derive, for a given SUSY spectrum, the bare parameters (according to Eqs. fl64|) and fl66l) ) 
that correctly reproduce physical quantities. 

After that, the effective couplings among quarks and leptons (in the mass eigenstates) 
with the Higgses are found performing the unitary transformations of Eq. fIBSl) on Eq. 

The procedure for the numerical calculation is summarized in Appendix O In the fol- 
lowing, we derive the analytical expressions for e{j (/ = u, d, e) and Uf{f = ul, ur, di, dR, 
CliGr) using the mass insertion (MI) approximation to explicitly show their parameter 
dependences on the flavor structures of the MSSM. A comparison between the numerical 
results obtained by means of the MI approximation and the full computation shows that 
the two approaches are fully consistent provided the MI parameters are small enough, as 
expected. 

The calculation of the loop induced non-holomorphic Yukawa couplings requires to 
specify the SUSY spectrum, which depends on the hare Yukawa couplings through the 
left-right mixing terms in the sfermion mass matrices. Therefore, it is useful to perform 
the calculation using the MI parameters defined in the SCKM basis for hare Yukawa 
couplings {hare SCKM basis). Clearly, physical observables do not depend on the choice 
of the basis and we express the MI parameters and the effective vertices for squarks in 
the hare SCKM basis, just for convenience. 

With an appropriate redefinition of the quark fields and the hare CKM matrix, the 
unitary matrices can take the following form, 

Ud,, = e-''^Va,, U,, = V,„ = e-^^~"K„ f/n, = K„ (67) 

where the diagonal components of Vf{f = dR,dL,UR,UL) are real. In the following, we 
refer to this basis as the hare SCKM basis. The formulae given in Section [3] are valid 
for the hare MI parameters defined in this basis. The MI parameters defined in the hare 
SCKM basis can be related to those defined in the physical SCKM basis by the following 
unitary transformation 

iSicx)':^'' = (ufAS^xxr-'^U}] , (68) 

V / ij 
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where X = L,R, f = d,u and i,j are flavor indeces. In practice, their difference can be 
safely neglected as long as the off-diagonal components of the unitary matrices U and 
their phase factors are small enough. 

In the SM, the quark and lepton mass matrices show an highly hierarchical structure 
and only the Yukawa couplings of the third generations induce large radiative correc- 
tions. To pick up them, we define a projection operator A = diag(0, 0, 1). Under the 
approximations mentioned above, the corrections for the quark and lepton masses are 
parameterized as 
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~~6 



^Rji^P 



i^RR)ij^ej 



LL)ijy 



(69) 



where the loop factors are given in Appendix [B] and a summation over the index k is 
assumed in Eq. fl69|) . Assuming equal SUSY masses, it turns out that eLR = = = 
= ^'lr = 4 = ^'r = ^'g = sign(/iM3) x a^/STT, ey = e^y = -sign(/iAt) x y1lZ2'K'^ and 
e'y = e'j^y = — sign(/iy4{,) x yl/?,2'K'^ for quarks. For leptons, = ei = — sign(/iM2) x 
80:2/1671, = sign(/iMi) x ay/Svr and = 0. 

Next, we derive the approximate expressions for the unitary matrices f/^^ and Ud]^ 
that diagonalize the down-type quark mass matrix. As stated before, it is convenient to 
choose the basis where 



Ud 



dL 



dL ' 



(70) 



where the phase rotation e ^^'^ has been introduced to make the quark masses real. From 
Eq. f lMl) . Udj^, Ud^ and e~*^'* are determined by the following matching condition. 



Expanding V^^^ 



Vdj^i^d + 5md)Vl = rhd 
we find 



mdi - md] 



mdi - rrid^ 



{nidi + {S'rn\)ii){5md)ij + {5m\)ij{mdj + {5md 
fUdi + {Smd)ii){Sml)ij + {6md)ij(jndj + {5m\)j 



JO) 



(71) 
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and 



where 



mdic'''^ = rridiil + eit/3 + efhD , (73) 



rrid. 



(2) 2 _ _JMdyMhi)^ , _ . 



In the expression of ^f^tj^-, we have systematically neglected subdominant terms of order 
0{\{5md)a\'^ /ml) because always negligible. Moreover, we have also neglected the second- 
order terms in the expansion of the unitary matrices V^^ and V^^; this approximation is 
valid as long as the MI parameters are significantly lesser than one. 

Notice that Eq. (!73|) determines the phase factor and relates hare and physical 
quark masses in the following way, 

"^di = (2) 01 ' = + + e^tj) • (75) 

\l + e,t0 + ert}\ 

Using the parameterization of Eq. fl69l) . we derive 

,r.r X (^R3itp ,^d \ "^di ZYi^^*tp—* nidi ZLjeLlitp ^ d \ (^a\ 
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9(1 + est/j) 3(1 + e-itp) md^ 

ey(eLj3 — eizi + exyj) ^ 



3(1 + e^tp) y LL In. K ) 



for i = 1, 2 and 



e3 = 6,3 + ey-5^^(\/^A5l,),., 
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The parameters zy,zr and zl appearing in the above expressions are defined in Ap- 
pendix [HI Although terms proportional to these parameters are suppressed by the light 
quark masses (and thus negligible in many cases), they are still relevant for the calculation 
of the (C)EDMs when these arise from J^l. 

The components of the bare CKM matrix can be expressed as a combination of physical 
CKM matrix elements and mass insertion parameters. In Eq. f l66l) . the contribution from 
Uu^ is negligible due to the suppression of t^^ as will be seen later. Thus, the components 
of the bare CKM matrix are given as 

1 + eitp 3(1 + e^tp) 



where 63 = es — ey and 



i'. = T^fl + #), (80) 



1 + eltp V 



where i ^ j with i,j = 1,2. Although pi ^ 0{etp) is comparable to other terms in 
Eq. fl79|) . we set hereafter pi = 0, for simplicity. In terms of the physical CKM matrix, 
Eqs. fITBl) and (1771) can be rewritten as 

^^^'^^^^^ - -W^)^ ''^'' ^ ITi^^^^ " ^3(l + e^t,)^^«^^='^' 

3(1 + est,) - ^T^^'^ + ^3(l + 6*t,)^^^^^^^' 

where and zr are defined in Appendix [B] and 

= + — (t^iifl^t^LL)^ 5-('^LL^>/)n, 

9(l + e3t^)m,^^ «^ 3(l + e3t^)m,^^ ' ^ ^ 



for i = 1, 2 and 



.=1,2 3(l + e^t,; 

^^^2^ 9(1 + 63^/3)^ ^« 3(l + e3t;3)^ ^ ^ 
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In the above expression for e] , we have omitted terms proportional to ev^La and 6^*6^3^ 
since they vanish for degenerate squarks and vanishing flavor-conserving CP phases. 

Turning to the up-type quarks, a similar calculation can be implemented and the re- 
lated expressions are obtained by the corresponding expressions for the down-type quarks 
by exchanging V {tp) with v'^ (^/3^)- We define unitary matrices as 

f/«H = e-*'"K„ f/., = K„ (84) 

and we match the hare Yukawa couplings with physical masses. In contrast to the case of 
down-type quarks, the corrections for up-type quark masses are suppressed by tan /3, 

= arg(l + e:t^i). 



(85) 



where 



9i ^ 



1,2) 



i=l,2 



LL)ii 1 



(86) 



and 63 = eg — Cy. In the above expression, while terms proportional to can be 
safely neglected, those proportional to t^^ should be retained. In fact, the e- parameter in 
Eq. (185!) is related to the vertex corrections for the up quark couplings with Higgs, leading 
to tan /3-enhanced corrections, as we will discuss in the next subsection. Concerning the 
unitary matrices, by means of the parameterization of Eq. (!69|) . we find 
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" mt 3 ^ 

where terms proportional to are safely neglected and z'y^z'^^ and 2;^ are defined in 
Appendix [Bl 

The same formalism used for the quarks applies to the leptons, too. In particular, the 
unitary matrices that diagonalize the lepton mass matrix can be defined as 



U, 



(88) 



and the phase factors, which make the lepton masses real, and the corrections to the 
lepton masses can be parameterized as 



rriei 
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arg(l + e^tp 



(89) 
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With the parameters in Eq. (!69l) . we find 



3(1 + e^tp) 
3{l + ef^ 



m^h^^UT^^i^RRh, (90) 



and 



= 4^+^^^iS%RA5l,h (^ = 1,2), 



and = e^^. 



A. 2 Effective Higgs Interactions 

We are now ready to derive the effective interactions of matter fields at the BLO. In this 
subsection, we consider the interactions of fermions with Higgs fields. The effective cou- 
plings in the physical CKM basis are obtained by performing the unitary transformations 
of Eq. (jHSD to the Higgs couplings with quarks in the bare SCKM basis, which is described 
by Eq. ([62]). 

For charged Higgs interactions, we obtain 

i-eff - —E, */3 I ) URidijH + — = tp I 1 ULidRjH 

+^^tJc^^) uueR,H+ + h.c., (92) 

where 

^^"^^^^^ = ^(U'^n^-^L^w + tl) - yvt% (93) 

''''ui 

^^Sh = -^iVU,,m,UlUl + t/) - V,,t/ , (94) 

''''dj 

In the above expressions, the non-holomorphic Yukawa terms e'fp^ are expressed by the 
Yukawa couplings and the unitary matrices through Eq. (16^ . then terms proportional 
to tl or t^'^ in Eqs. (EMM]) come from the one-loop vertex corrections for charged Higgs 
couplings while the remaining terms are generated by the mass term corrections. Notice 
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that, in the hmit oi Uf — 1 {f — dL,dR,UL,UR,eL,eR) and m/ = m/ (/ = u,d,e), the 
vertex corrections are vanishing, as it should be. This means that the vertex corrections 
are encoded in the difference between bare and physical Yulcawa couphngs and in the 
deviation of unitary matrices from the identity matrix. Especially, we need to keep tan (3- 
suppressed terms in Uu^ , and m„ for C^^ . 

Neutral Higgs interactions with quarks and leptons are obtained in a similar way by 



where 



= Z^iUu.^uUlha^* + 5,,sphT , (97) 



m. 



iCa:h = --^iUa,'rn,Ul),ja^ + S,,cpbr , (98) 

(^2")^^- = -:^(Ue,'^eUlJ^Ja'^ + S,,C^br , (99) 

with a" = 0^62* ~ Sfjb^. The Higgs fields and iJ^ are expanded as hnear combinations 
of mass eigenstates Ha — h, H, A, and the coefficients 6" and 63 are defined as 

Hl^^{vcp + blHa), Hl = ^{vsp + blHa). (100) 

If there is no CP violation in the Higgs sector, b1 = {—Sa, Ca, ise}, b^ = {cq,, Sq, icp}, 
= {Ca-i3, Sa-pi — i} and then we can identify h and H as the CP-even Higgs bosons 
and A as a the CP-even Higgs boson. In the presence of CP violation, mass eigenstates 
can not be assigned with CP charges, and all of the above coefficients become complex 
numbers as 



bl = Ola + ISpO^a, 6^ = 02a + ICpO^a, (101) 

where the derivation of the matrix Oia can be found in Refs. [38,39]. 

To calculate the physical observables keeping their dependence on the flavor- and CP- 
violating parameters in the MSSM, it is useful to derive analytical approximate formulae 
for these vertices. Since the unitary matrices and tree-level Yukawa couplings are already 
derived by means of MI parameters and loop factors, we can easily derive the expressions 
for the desired vertices. Here, we present only the vertices that are relevant for the 
calculations of the flavored EDMs. 

In the calculation of quark (C)EDMs, flavor- violating quark couplings with the charged 
Higgs have a very important role, as shown in Ref. [12]. The corrections to the Higgs 
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couplings with quarks, induced by the non-holomorphic Yukawa couphngs, are given as 
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(103) 



for i = 1,2. Here terms proportional to z'x {X = R, L, Y) are neglected since they 
are suppressed by rriujint, thus, irrelevant. The last term in Eq. fll02p is significantly 
enhanced by the bottom Yukawa coupling (through the mass corrections for down-type 
quarks); the enhancement factor (mb/m^J and the appearance of an extra factor par- 
tially compensate the suppression arising from the loop factor and the MI parameters. 
As a result, these fiavor- violating terms can become comparable in size to the tree-level 
couplings and with a different phase. 

For light up-type quarks, the couplings with the charged Higgs are given by 
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(105) 



for i = 1,2. In the above derivation, we have assumed {S^lV)-^ ~ (5LL^)3j — ('^ll)3* 
and terms proportional to z'j^ {X = R, L, Y) have been neglected. The tan /3-enhanced 
corrections appear from the vertex corrections in contrast to the case of down-quarks, 
where the corrections come from the mass corrections, and especially the last term of 
Eq. ( 11051) is enhanced by the large top quark mass. 

For the neutral Higgs interactions, fiavor-violating interactions are expressed in a 
similar way. Moreover, fiavor-conserving CP-violating interactions appear at the BLO 
and they are important for the calculation of Higgs mediated CP-violating four-Fermi 
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interactions. In the limit of large tan/3 and niA, these couplings are expressed as 



( ef\l 



(2)2, r 1^.. ' V— 1'2), (106) 

(C,^J,, ^ (107) 

(CJ,, ^ -c2(l-e:t^), (108) 

(C^ik ^ z4(l + e:t^), (109) 

and {C^j)ii — iC^j)ii — spCp since the lightest Higgs boson behaves as the SM Higgs 
boson in this limit. As we can see from Eq. fll06p . CP- violating interactions between 
down-type quarks and the Higgses H and A can arise from i) diagonal corrections to the 
quark mass matrix (e~*^'), an effect which has been already discussed in the literature, 
or a) from the combination ((V^^AV^^^)jj), i.e. from pure flavor-violating effects to the 
quark mass matrix. 

Although these last effects are suppressed by small flavor- violating parameters, they 
are highly enhanced by the heaviest quark Yukawas thus, they can provide potentially 
large effects. 



A. 3 Effective Fermion-Sfermion Interactions 

The effects of non-holomorphic corrections appear not only in the interactions with Higgs 
bosons but also in those with the Higgsino and gauginos. Here we discussed the effective 
vertices including the BLO tan /3-enhanced corrections. 

First of all, we need to specify the basis in which we work. In our convention, while 
quarks are expressed in the mass eigenstates, we adopt the hare SCKM basis for squarks 
and flavor eigenstates for gauginos and Higgsinos to keep the dependence of flavor and 
CP-violating parameters. Thus, the left-right mixing terms proportional to the yU term 
in the squark mass matrices become flavor blind (see Eqs. fll33p and (11361) ) while BLO 
corrections induce flavor-violating couplings even for neutral SUSY particles such as the 
gluino and neutralinos. If we use physical SCKM basis, these couplings with neutral 
particles become flavor blind but there appear non-trivial flavor violations in the left- 
right squark mixing terms. Obviously, physical observables do not depend on the choice 
of the squark basis. In Appendix O, we also provide the relevant vertices for all particles 
in mass eigenstates. 

First, gluino interactions are modified by the BLO corrections as 

- 4„t = {Cl)^^ TjltaqL, + {C'lh QR^LtaqR, + h.C, (110) 

where 

(Clh = -^9s{Ul)., (Ill) 
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with q — u,d. The flavor-mixing terms are generated by the mismatch between the 
bare and physical SCKM bases. It is worth mentioning that even the flavor- conserving 
couphngs can have CP- violating phases at the BLO due to the phases 9q in the unitary 
matrix t/^^. 

The charginos are the mass eigenstates of charged Higgsino and wino, which are mixed 
each other after the electroweak symmetry breaking. Including the BLO corrections, the 
vertices of these charged Higgsino and wino are given by 

+{CZ%J^uW^R^Lj + {Ci%je*^MR^Lj 

+ h.c. , (112) 

where 





)i3 


= ^2([/i,n.-, 




lij 






hj 






lij 






lij 






lij 






kj 




(c-r 








lij 


= fl'2%, 




lij 


= iVe)ij ■ 



The couplings of neutralinos, which are composed by bino, wino and neutral Higgsinos, 
are given by 

- ^fnt = i^Z )ijdLiWRdLj + {C^Jijdl^BRdLj + {C^j^)ijd*jiiBLdRj 
+(C<^ )ijd*RiH2RdLj + (C^ )ijd*uHl^dRj 

HcZ\^LwWj + ic!J^Jirl^L, + {cihu*R^LUR, 

+ icT)A^Wj + {Cl),fel^^eL, + {Cih~e*R^LeRj 
+(CeL )ij^*RiHmeLj + (C^j^ )ije*LiH^LeRj 

+iCZ\i^lW>Lj + {cihvl^RVLj + h.c. , (114) 
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where 



dL 



{C. 



HO 
H° 



(nilo 



M/0 



^92 



2 y'^UL^v 



V2gyY,,{Ul)v, 



V2g2 



V2gyYiMjij, 
V2gyY,Mnh , 

V2g2, 



V2gyYiJij , 



(115) 



where Y,^ = 1/6, Y^, = 1/3, = -2/3, F;, = -1/2 and Y,^ = 1. 

It is easy to see the parameter dependence of these effective couphngs just considering 
the approximate expressions for the unitary matrices. For example, flavor- violating gluino 
interactions with down-type quarks can be written as 



{CDs, c V2g, 



3(1 + est 



(C|j3^ -^9s e 



id,. 



l + est 



/3 



rfik 3(1 + eltp) 



RR)3i 



, (116) 



3(1 + eltp) ^^^^^'^ 



mdi Zy*eyt(3 



. ^di ZL*eL3it0 , ^ , 



rUb 3(1 + estf)) nib 3(1 + 63^/3) 



;ii7) 



for i = 1,2, and the flavor conserving couphngs of the right-handed down- type quarks 
have CP phases as (C^^)ii = —V^gs e^^'^t, in our convention. 
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To discuss the down-type quark (C)EDMs, it may be convenient to take the basis of 
left-handed up squarks same as that of left-handed down squarks. In that case, effective 
couplings and C^^ should be replaced by C]^° and Cj^", respectively, and also the 
left-right mixing mass term of up squarks is given by 

- £ = {ml^)ij UR-ULj + h.c. , (118) 

where 

{m\L)i3 - -mui{Au + fi* cot (3) {U^^VUl)i, . (119) 

Here, the A terms are assumed to be proportional to the Yukawa couplings as {au)ij — 
A^y^^ in the hare SCKM basis. 



B Loop Factors 

In Appendix \^ non-holomorphic Yukawa corrections are expressed by the combination 
of mass insertion parameters and loop factors. Here, we present the expression for the 
loop factors. 

For down quarks, the loop factors are given by 



a. 



^*M; h{\M^\\m 



_ _ _ _ _ Tfl ~ 

:/x*M*[2/3(|Ml|^mi.,ml.) + 3/3(|Ml|^|/i|^m^- ,) + 6/3(|Ml|^|/i|^ml.) 



%*M*J3(|M2|M/ir,ml,) 



*7i/r* <->T/i7i^|2 ^2 ^2 \ I o T l \ 'h/T \'i l,.|2 ^2 \ I ar l\ 'h/r \2 l,.|2 ^2 



(-Li 



727r 



2 2 

dm ' d^i 



dhi' 



dm 



q I ' dm' dLi' dLj' IQtt ^ '^'^ " 



dhi ' rfij ' 



(Rij — 



^LRij — 



— u*M*m|/4(|M3|2,m 

3^^ 3 d 4VI 31 , 

ay 



222 

dLi am ajji 



22 2 ' 

affi ' dnj ■ 



a 



^/i*M>X- J5(|M3|^m 
37r " 



ey = 



2 2 2 2 2-1 

'^";t ,"017 ,"0^7 ,"017 , 



/z*M>X- /5(|Mir,m| ,m? ,m5- ^,m? ) 



f2 



CLyj — 



0:2 vai 
IGvr 



1 + -2 ) m? /4(|/i|2,m?^,m?^,m^^J . 



(120) 



The functions /3, J4 and J5 are defined in Appendix [El 



43 



For the up quarks, the loop factors are 



2 



ay 



-Lij 



1447r 

s 

37r 



f^^*M* 4J3(|Mi|2,mi,,,mL) + 3/3(|M2p,|/.p,mL)-12J3(|M2p,|/xp,m^ 

3q;2 



2 



'^s * ht* 2 t n n/T 12 2 2 2 \ ^"2 *ii/f* 27/171^121 |2 2 2 



IGtt 



+ 



ay 



^*M*mi 4/4(|Ml|^m^^^,m^^^,m|^^) + 3/3(|M2|M/x|^m^^^,m^^^) 



a 



'Rij 



ay 



367r 



ay 
36n 



2 _2 _2 



0:2 
IGtt 



+ /3(|/.|^ m?^,m?J, 
(1 + m? hi\i^\'^,ml^,ml^,mlj . 



(121) 



Concerning the expression for the off-diagonal components of the unitary matrix that 
diagonalizes the quark mass matrices, it is useful to introduce some combinations of the 
e parameters. For down quarks, we define 





1 + £3^/3 


ZRi = 


1 + Ut/s 1 

1 + €3^/3 


ZLi = 


l + Qt^ 

1 + estf, 




1 + £3^/3 


ZLi = 


1 + eit/3 1 
1 + £3^/3 



eL3i 1 + £3^/3 1 + £3^/3 

^Ri3 1 + £3^/3 



(122) 



and for up quarks we define 



~ 1 + 



~ 1 + 



--Ri3 
^/ * ' 
-i?3i 



^Li — 



1 + 



--Li3 
c' * 



(123) 
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The loop factors for charged leptons are calculated 



as 



4 = 



^/i*M* /3(|M2p, |/ip,m^^J + ^^,*Ml\2h{\Ml 
+J3(|Mi|M/i|2,m?^J-2/3(|Mi 



_^|2 |,.|2 



I |2 2 N 



lovr 

tty „ r 



|2 I |2 2 2 \ 

!2| ,mg^^,mg^^O 



+ ^/^*Mrm? [2/4(|Ml|^m^^^,m^-^^,m^^^.) + /4(|Ml|^ 
= ^^*M>^ [/4(|Ml|^m^^,^,m^„,,m^„J-/4(|Ml|M/i|^ 



ay 

where the functi 



I 1 2 2 

■11 , m 

2 2 2 \ r /I ji^ |2 I |2 2 2 \ 

,m= , — J4 Mi , u ,m= , 

ELi' efli' eflj' -'-I ' lA^I ? e^ii' eflj' 

ions I3, I4 and Is are defined in Appendix [El 



(124) 



C Effective Vertices for Mass Eigenstates 

In this Appendix, we present the numerical procedure to derive the effective vertices for 
the scalar mass eigenstates, which are used in our numerical calculation of Section HI 

First, we need to derive the non-holomorphic Yukawa corrections for a given set of 
tree-level Yukawa couplings and soft SUSY-breaking parameters. To this end, we have to 
diagonalize the sfermion mass matrices by means of the unitary matrix Uj^ as 



Uj m) U) 



m ; 



if = q,u,d,l,e) . 



(125) 



Using the unitary matrices and tree-level Yukawa couplings yu = yu^^ Vd^Ud ^'^^ Ue^Ve 
the non-holomorphic corrections defined in Eq. ( l62l) are described as follow 



(e-v^). 



"V*M3*(t/i),fc(f/,y,[/j)fcKf/.-)^i/3(|M3|2,m^,,m2 



11 



327r2 

3^2 * 



3l|/xp,m?^,m| 



167r 

ay 

'487r 
ay 

ay 
36n 



/i*M* (f/i),fc(f/,yJfc,/3 (|Ml|^ \^l\\ml^ 



fx*M* {UlUiUuVjJhkimi.h (\M,\\ml,,m, 



(126) 
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03 



3^2 



167r 
ay 



7^ 



^^2* (y,f/i).fc(f/,-).,/3 (|M2|^ |/i|^mJ-^ 



487r 



(127) 



4 = ^l^*M;iy,UlUmkjh{\M2\^\fi\^ml) 



(128) 



The loop function I^^x, y, z) is defined in the Appendix [El 

Mass eigenstates for quarks and leptons are determined by Eq. flMl) and we can nu- 
merically derive the unitary matrices Uf {f = dL, dn, ul, ur, cl, en) which relate the mass 
eigenstates and fiavor eigenstates. In the procedure, tree-level Yukawa couplings must be 
iteratively adjusted to reproduce experimentally-observed values of the CKM matrix and 
mass eigenvalues. After the achievement of the adjustment with the desired accuracy, 
tree-level Yukawa couplings are determined, which fix the SUSY spectrum, and finally 
unitary matrices Uf are obtained. 

Although the effective vertices among quarks/leptons and SUSY particles have already 
been derived in terms of Uf and yf in Appendix [XJ we derive here these vertices also in 
the mass eigenstates, for numerical convenience. 

Firstly, unitary matrices that diagonalize the mass matrices of SUSY particles must 
be defined. For the sfermions, the unitary matrices f/j are given as 



Uf 



m 



RL 



m 
m 



2t 
RL 
2 

RR 



U) = rn)., 



(129) 



with f = d,u,e and mass eigenstates // and fiavor eigenstates in the bare SCKM basis 
fi = {fi, Ir) are related by 



fi = {Uj)jJ, (J 



,6, i 



.6) 



(130) 
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In our convention, for up quarks 



{ml,) 



{Vm~V')ij + ( ^ -^sin^Ow ) m| cos 2/3 5^ , 



[m 



RRjij 



[m 



RLjij 



and for down quarks 



(m 



RRjij 



+ 2 sin 9w cos 2/3 %, 



= (m?)y + (^~\ + \ si^^ ^w^) "^1 cos 2/3 , 



(m|)ij - - sin^ 6w m| cos 2/3 % , 

V cos /3 



V2 



in the bare SCKM basis. The mass matrix for charged sleptons can be read 



{m~)ij + + sin^ 9w^ rn% cos 2/3 

(m|)jj — sin^ 9w fn\ cos 2/3 5^ , 

V cos 13 I _ . * \ 



while that for left-handed sneutrinos is defined as 

1 

{m\j)ij = (mf + -m| cos 2/3 (5^^ , 

and it should be diagonalized by a 3-by-3 unitary matrix. 
The chargino mass matrix is diagonalized by 



O 



M2 \/2mw cos /3 
^ V V^fnw sin /3 /i ' 



Ci 



C2 



Charginos are mixed states of charged wino and Higgsino defined as 



XlL 




X2L . 




XlR 


^Or 


X2R 



H 



IL 



H. 



R 
2R 



and the charge conjugate fields are defined by V'L/ij — {'^^r/lT '^'^ ~ 
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u 



N 



The neutralino mass matrix is diagonalized by 

f Ml -mzswcp mzswsp \ 

M2 mzCwCjS -mzCwSf3 

-mzswcp mzcwcp — /i 

\ mzswsp -mzcwsp -/x J 



( 



N 



M 



N2 



\ 



Af4 



This unitary matrix is easily obtained by considering the square of both sides of Eq. fll43p . 
Neutrahnos are mixed states of bino, neutral wino and Higgsinos defined by 



(144) 



and xm = 

Using the above unitary matrices and the effective vertices in the flavor basis given in 
Appendix. [Aj which are expressed by tree-level Yukawa couplings and unitary matrices 
Uf, one can find the effective vertices in the mass eigenstates. 

The couplings of the gluino to the quarks are given by 



£L = rrgX [{^Ji^Pl + (rL)/,P«] q, + h.c. 



where 



k=l 

3 



(145) 

(146) 
(147) 



fc=i 



The couplings of the chargino to the quarks and leptons are given by 



+(rfJt,^>;x^ie,, + (rfj^^.z>;xZe/^, + (r^j,V~/x^i^^L, + h.c. 



(148) 



where 





= E [iOR)kimniCZ% 
1=1 


+ {OR)k2{Uu)li+3{C^^ )ij 




3 

1=1 






= E [ioi)ki{uMcZ% 

i=l 






3 

= J2(OR)k2{uMc!!^h, 





(149) 



i=l 
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for quarks and 



1=1 

3 

J2{OL)k2mn{C, 

i=l 
3 



'en Jij 5 



i=l 



(150) 



for leptons. The couplings in the flavor eigenstates are defined in Eq. (11131) . 
The couplings of the neutralino to the quarks are given by 



(151) 



where 



(rg)/, = E + iUnU{uMcl\ + (f/^).3(f/j)/.+3(C;: 



i=l 
3 



(ri)/, = E Yu%U{u^)uUcl\., + (f/^).3(t/j)/.(<).. 

3 

(r£)/, = E + {UM),,mn{cZ\ + ([/^)M(f/.)/.+3(cf;) 

3 

(r£)/, = E [(f^^)'^i(^^)/^+3(c4)., + {w^)m{Uu)h{c, 



'UR Jij 



(152) 



i=l 

for quarks and 

3 



i=l 

i=l 
3 

(r£)/, = [E(^^)'^i(f^^)^^(<)^^ + (^^)'^'2(f^'>)^^(C', 

i=l 

for leptons. The couplings in the flavor eigenstates are defined in Eq. (11151) . 



(153) 
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D (C)EDMs Formulae for Mass Eigenstates 



Here, we present the full expressions for the quark (C)EDMs and lepton EDMs in the 
mass eigenstates basis using, in particular, the effective vertices derived in the previous 
section. 

The gluino-mediated contributions read 
{d^^^ dl}, - E [Ms (ri),:(rl),] {/f (x,) , /f (x,)} , (154) 

where Xj — |M3p/m|^.. The chargino-mediated contributions read 

di),^ = -3^5 E E H^^:'xtn] {jp (4) . fP . (155) 

j=i k=i '^j 

where a;^ = MQ^/m\-. The neutralino-mediated contributions read 

1 ^ 



rrir 

j=l k=l dj 



where x'^ — M'^^/m'^-.. Finally, the Higgs-mediated contributions have the following 



expression. 



[dd. dl) 



2 3 ^2 
92 ^d. 



647r2 ^ m?,, 



(d.Ad..] (^.) , ^ (^.)} , (157) 



where Xj = ml. / . 

Passing to the up quark (C)EDMs, the gluino-mediated contribution is given by 

{du. = E ;^Im [M3 (r|J,/(rU,,] {/f (X,) , /f (x,)} , (158) 

j=i 

where = |M3p/m|^-. The chargino-mediated contribution reads 



6 2 



{rf... < } J. = - 3^ E E [(r£)f;(r£)J.] {/f (4) . /r (4) } . (159) 

j=l k=l dj 



where Xj — M^j^/rn?^.. The neutralino-mediated contribution reads 



{<• = -iEE^i- [(r£)r(r5:)j] {tf' (4) • /r (4)} • (m 



j=i fe=i 
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where Xj = Mf^^/ml^-. Finally, the Higgs- mediated contribution is given by 



92 



2 — ml 



H± 647^2 ^2 ^2 



[C^UCth*] {f'S (^.) , f? (^.)} ,(161) 



where Xj = ml/M]j±. 

The lepton EDMs receive contributions by the only chargino and neutralino sectors. 
The chargino-mediated contribution reads 



1 A A M, 



j=i k=i "J 



where Xj = Ml;^/m\-. The neutralino-mediated contribution is given by 



^" 327r2 



1 ^^^i^ 



j=l k=l ^3 



(rg)^ (rS)J. (4) 



where = Mj^j^/ml-. 



(162) 



(163) 



E Loop functions for EDM Formulae 

The loop functions that appear in Sections [2] and [3] are defined as 

1 — + 2x log X 



/r(-) 



(1-X)3 

3 — 4x + + 2 log X 

(1-X)3 ' 

— 1 — 4x + — 2x(x + 2) logx 
(1 -x)4 



-5 + 4x + x^ - 2(1 + 2x) log X 

(l-x)4 ' 
1 + 9x — 9x^ — x'^ + 6x(x + 1) logx 

(1-X)5 ' 

2(3 - 3x2 ^(1^4^^ 2,2) iQg^) 

(1 -X)5 ' 

-3 - 44x + 36x2 ^ ^2x3 - x^ - 12x(3x + 2) logx 
3(1 -x)6 

-10 - 9x + 18x2 + x^ - 3(1 + 6x + 3x2) j^g^ 
^ 3(1 -x)6 ' 



(164) 
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and 



X - y 



(i,^ = 0,l). 



In the limit of degenerate arguments, the functions are given as 

/o^ni) = <i^, i?;, -i^} (z = 0,l,2,3), 

fi\^) = i-^, i:, -i^, i^} = 0,1, 2, 3), 



111 1 

3' ~ 6' 10' ~ 15 

2 1 111 
~3' 6' ~ 15' 30 J 

r 1 ^1 

\ 6' 15/ 
#(1,1) = -i| (j=0,l). 



The following combinations of functions are also used 



-3/^\^)+/r(^), 



(165) 



(166) 



(167) 
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F Loop functions for Yukawa coupling corrections 

The loop functions used in the derivation of non-holomorphic Yukawa couphngs are given 

by 

/ aloga blogb clogc 
73(0, t), c) = -2 + — + 



[c — a){a — b) {a — b){b — c) {b — c){c — a)J ' 

IJa, b, c, d) — —6 (7 oga ^ ^ cydic permutations ] , 

\[a — b)[a — c) [a — a) ) 

hia, b, c, d, e) = 12 ( — r-^ — + 4 cydic permutations ) . (168) 

\{a — b){a — c) [a — d)[a — e) ) 

These functions are normahzed as /3(1, 1, 1) = /4(1, 1, 1, 1) = -^5(1, 1, 1, 1, 1) = 1- 
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